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陈 


肝 内 影响 最 为 广泛 的 自然 灾害 之 一 ,其 所 导致 的 土壤 沙漠 化 、 萝 漠 化 和 熙 碱 


化 给 生态 环境 造成 不 可 闭 的 危害 。 通 过 对 MODIS 数据 进行 投影 转换 、 去 云 等 预 处 理 的 基础 上 , 利 
用 地 形 校正 对 TVDI 模型 进行 改进 ,构建 了 改进 型 的 温度 植被 干旱 指数 (m7TVD1) 用 于 新 疆 干 时 区 旱 


情 监 测 。 利 用 土壤 实测 数据 对 mTVDI 及 传统 的 


TVDI 模型 进行 对 比 验证 。 研 究 结果 表 明 :(1) Fl 


用 EVI 与 校正 后 的 LST 构建 的 mTVDI, 对 干旱 区 旱情 的 敏感 度 最 高 ,与 实测 土壤 水 分 数据 的 相关 


KE R^ 为 0.74。(2) 从 空间 上 看 ,新 疆 2015 年 旱 


心 ,旱情 状况 由 严重 逐步 向 周围 山 


情 分 布 以 塔里木 盆地 和 准噶尔 盆地 为 两 个 干旱 中 
区 说 减 至 温润 状态 。 从 时 间 上 看 ,新 疆 6 月 7 月 和 8 


月 旱情 最 


为 严重 。(3) 研究 利用 TRMM 降水 数据 对 基于 mTVDI, 反 演 的 新 疆 旱 情 时 空 分 布 特征 进行 对 比分 
析 , 结 果 表明 二 者 所 表现 出 的 旱情 时 空 分 布 较为 一 致 ,不 同时 间 段 内 的 降水 量 与 m7TVDIs 之 间 具 有 
一 定 的 相关 性 , 且 均 通过 了 已 <0.01 显著 性 检验 。 综 上 ,基于 TYD7 所 提出 的 mTYD 大 能 够 有 效 开 展 


新 疆 干 旱 区 旱情 监测 , 且 精 度 较 高 ,从 而 为 今后 定量 化 开展 大 区 域 尺度 旱情 监测 研究 提供 参考 。 
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干旱 是 全 球 范 围 内 最 常见 .最 复杂 的 自然 灾害 
之 一 ,具有 影响 范围 大 .持续 时 间 长 .破坏 性 大 等 特 
点 ,不 仅 造成 巨大 的 经 济 社会 损失 ,同时 对 生态 环境 
产生 重大 威胁 。 近 年 来 , 受 全 球 气 候 变 化 的 影响 , 干 
早 区 旱情 状况 更 加 复杂 ,因此 定量 化 开展 大 区 域 尺 
度 旱情 监测 研究 对 于 把 握 全 球 干旱 区 气候 变化 规 
E ,保障 生态 安全 、 预 防 旱灾 的 发 生发 展 具有 重要 的 
现实 意义 。 传 统 旱情 监测 以 现场 监测 方式 为 
主 ,包括 : 探 针 测量 .时 域 反 射 测量 等 诸多 方法 ,但 是 
这 些 方法 消耗 巨大 的 人 力 物 力 且 效率 低 , 尤 其 是 山 
区 和 偏远 地 区 采样 困难 ,使 得 大 范围 监测 难以 实现 。 
遥感 技术 因 其 具有 监测 范围 广 、 光 谱 特 征 丰富 等 优 
势 , 有 效 弥补 了 地 面 站 点 监测 的 不 足 ,已 成 为 目前 旱 
情 监 测 的 重要 手段 "” 。 当 前 以 TVDI( Temperature 
Vegetation Dryness Index) 为 代表 的 植被 指数 遥感 旱 
情 监测 方法 ,通过 利用 植被 冠 层 温度 和 植被 指数 的 
变化 特征 ,构造 出 地 表 温 度 一 植被 指数 的 特征 空间 ， 
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耦合 温度 植被 干旱 指数 ,能够 对 不 同 下 垫 面 实现 大 
范围 长 时 序 监测 ,已 在 旱情 监测 方面 得 到 广泛 应 
Jj -*', SCHIRMBECK 等 551 比较 了 两 种 TVDI 模 
型 用 于 巴西 大 豆 种 植 区 的 旱情 监测 ,结果 表明 基于 
作物 参数 化 构建 的 TVDI 模型 能 够 确定 作物 旱情 发 
生 的 时 间 和 频率 。ZHU Wenbin 等 “提出 了 一 种 改 
进 型 TVD7, 采 用 连续 干旱 监测 的 方法 突破 了 传统 依 
靠 空 间 分 辩 率 主导 方法 的 限制 。YANG Ruowen 等 
人 对 TVDI 干 边 拟 合 方程 进行 改进 ,利用 二 次 方程 
来 表示 干 边 拟 合 方程 ,对 2010 一 2012 年 云南 旱情 状 
况 进行 监测 呈 ] 。 国 内 众多 学 者 利用 TVD 模型 对 我 
国 旱 情 状况 也 做 了 不 少 工作 。 刘 艳 伟 等 利用 Land- 
sat 数据 构建 了 由 归 一 化 植被 指数 和 地 表 温 度数 据 
建立 的 二 维特 征 空间 , 拟 合 了 干 , 湿 边 方程 ,进而 计 
算 了 7YD1, 对 朝阳 县 旱情 状况 进行 了 分 析 和 评 
Mr? 。 葛 少 青 等 利用 三 种 干旱 监测 方法 对 干旱 区 
沼泽 湿地 土壤 水 分 进行 监测 ,指出 三 种 干旱 指数 与 
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土壤 实测 数据 之 间 存 在 明显 的 负 相 关 关系 “。 然 
而 这 些 工作 都 侧重 于 小 区 域 短 时 效 的 研究 ,针对 干 
早 区 大 区 域 尺 度 植被 生长 季 早 情 状 况 的 研究 较为 少 
The 

鉴于 此 ,本 文 以 MODIS( Moderate Resolution Im- 
aging Spectroradiometer ) 光谱 数据 为 数据 源 ,对 传统 
的 TVDI 旱情 监测 模型 进行 改进 ,将 改进 过 后 的 温 
度 植被 干旱 指数 (modified Temperature Vegetation 
Dryness Index , mTVDI) 用 于 我 国 新 疆 干 旱 区 植被 生 
长 季 旱 情 监 测 研究 ,以 期 为 今后 进一步 开展 干旱 区 
旱情 监测 提供 新 思路 。 


1 研究 区 概况 及 数据 源 


1.1 研究 区 概况 

新 疆 干 旱 区 介 于 73°20’ ~ 96°25°E, 34°15! ~ 
49°10'N 之 间 , 深 居 亚 欧 大 陆 中 心 ,温带 大 陆 性 气 
候 。 由 于 特殊 的 地 理 位 置 , 主要 依靠 大 西洋 湿润 气 
团 带 来 降水 , 南 性 年 降雨 量 20 ~ 100 mm , 北 性 年 降 
雨量 100 ~500 mm ,气温 年 .日 较 差 大 ,光照 充足 AE 
发 强烈 ,是 全 球 典型 的 干旱 半 干 旱 气 候 区 。 新 疆 地 
形 轮廓 为 “三 山 夹 两 盆 ” ,北部 阿尔 泰山 脉 ,南部 昆 
仑 山脉 ,天 山 山 脉 横贯 中 部 ,塔里木 盆地 和 准 格 尔 盆 
地 夹 在 其 中 。 全 国 第 一 第 二 大 沙漠 塔克拉玛干 沙漠 
和 古 尔 班 通 古 特 沙漠 分 布 于 此 ,自然 环境 复杂 多 样 。 
根据 新 疆 植被 物候 时 空 变化 特征 及 新 疆 2015 年 
植被 生长 实际 情况 ,2015 新 疆 植被 生长 季 为 4 ~ 10 
月 份 即 第 97 天 ~ 第 305 天 。 
1.2 数据 源 与 预 处 理 
1.2.1 适 感 数据 ”遥感 数据 主要 采用 MODIS 数据 
产品 .SRTM 高 程 数据 以 及 TRMM 降水 数据 。MO- 
DIS 数据 产品 具有 禾 盖 范围 广 、 时 间 分 辨 率 高 等 特 
点 ,能 够 满足 时 序 研 究 的 需要 。 本 文选 取 的 MODIS 
数据 产品 包括 16 d 合成 植被 指数 ( Vegetation Index, 
VI) 产品 MODIS13A2 和 8 d 合成 地 表 温 度 ( Land 


Surface Temperature , LST) 产品 MODIS11A2?" , # 
tH, MODISI3A2 包括 NDVI ( Normalized Difference 
Vegetation Index ) 和 EVI ( Enhanced Vegetation Index) 
两 种 数据 产品 。MODIS11A2 产品 为 8 d 合成 产品 ， 
需要 对 其 进行 进一步 处 理 , 合 成 为 16 d 产品 ,以 确 
保 与 植被 指数 产品 时 序 相 一 致 。 数 字 高 程 模型 数据 
( Digital Elevation Model, DEM ) 选取 SRTM ( Shuttle 
Radar Topography Mission ) fj fe Zi s , BUA 72 s ya Fal 
广 、 精 度 高 应 用 广泛 等 特点 ,经 重 采 样 后 分 辩 率 为 
1 km? 。 降 水 数据 采用 TRMM( Tropical Rainfall 
Measuring Mission ) 卫星 数据 TRMM3B43, 是 由 
TRMM3B42 数据 NOAA ( National Oceanic and At- 
mospheric Administration ) 气候 预测 中 心气 修 异 常 系 
统 的 全 球 格 点 雨量 测量 器 资料 和 全 球 降 水 气候 中 心 
的 全 球 降水 资料 合成 的 全 球 网 格 降水 数据 ,单位 为 
整 月 平均 每 小 时 降水 量 (mm h ”) ,空间 分 辩 率 为 
0.25。x0.25。, 经 重 采 样 后 分 辨 率 为 1 km ” 。 所 需 
ERAGE EAM fa ILE 1 所 示 。 

1.2.2 实测 数据 实测 土壤 含水 量 数据 来 源 于 
2015 年 6 月 野外 实测 采样 工作 ,共计 50 个 样 点 , 采 
样 点 位 置 是 根据 新 疆 不 同类 型 植被 分 布 情况 统计 抽 
样 确定 ,具体 位 置 如 图 1。 采 样 时 间 要 求 与 卫星 过 
境 时 间 一 致 ,天 气 要 求 晴 度 且 风力 较 小 。 采 用 土 外 
法 采集 0 ~ 10 em 土 层 样本 并 取 三 次 重复 ,用 铝 盒 封 
装 。 将 土壤 样本 带 回 实验 室 用 人 烘 干 法 人 烘 干 称 重 得 到 
土壤 含水 量 数据 。 具 体 公 式 如 下 : 
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W--a ¥ 100 (1) 
Wy 
式 中 :w 为 土壤 含水 量 (% ) ,wi 为 湿 土 质量 (g) , wy 
为 烘 干 后 土 的 质量 (g) 。 


1.2.3 数据 预 处 理 MODIS 数据 产品 需要 利用 
MRT( MODIS Reproject Tool) 工具 对 其 进行 几何 校 
正 .投影 转换 .格式 转换 数据 抽取 等 前 期 处 理 ,再 利 
用 AreccTS 对 研究 区 进行 裁剪 .波段 运算 合成 等 工 


R1 遥感 数据 详细 参数 信息 
Tab.1 Detailed parameter information of remote sensing data 

Co ee RN My "me 

数据 类 型 数据 名 称 数据 选取 时 间 ox 空间 分 辨 率 数据 来 源 
地 表 温 度数 据 MODIS11A2 2015 年 第 97 ~305 d 8 d ] km x1 km 美国 国家 航空 航天 局 https :// modis. gsfc. nasa. gov 
植被 指数 数据 MODISI3A2 2015 年 第 97 ~305 d 16 d 1 km x1 km 美国 国家 航空 航天 局 https ;// modis. gsfc. nasa. gov 
高 程 数据 SRTM 30 mx30 m ”美国 太空 总 署 美 国 国防 部 国家 测绘 局 ,德国 与 意 大 

利 航 天 机 构 http ://srtm. csi. cgiar. org 

降水 数据 TRMM3B43 2015 年 第 97 ~305 d 1 h 0.25? x0.25° 美国 气候 预测 中 心 https ;// mirador. gsfc. nasa. gov 
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Fig.1 Sketch Map of Study Area 


作 。SRTM 数据 需要 利用 ArcGIS 进行 裁剪 投影 
换 和 重 采 样 工 作 。TRMM 数据 需要 进行 裁剪 、 重 采 
样 和 波段 运算 。 


2 研究 方法 


2.1 植被 指数 

为 精确 反映 植被 信息 ,建立 植被 指数 与 温度 之 
间 的 关系 ,本 文选 取 归 一 化 植被 指数 NDVI 和 增强 
型 植被 指数 EVI 用 来 表示 植被 生长 状况 信息 。 
NDV1, 即 归 一 化 植被 指数 ,用 以 表现 植物 生长 状态 、 
植被 空间 分 布 密度 及 植被 冠 层 的 背景 影响 等 ,在 利 
用 遥感 图 像 进行 植被 研究 中 得 到 广泛 应 用 。EV, 增 
强 型 植被 指数 ,增加 了 有 冠 层 背景 调整 因子 ,使 得 EVI 
对 大 多 数 冠 层 背景 不 敏感 ,提高 植被 信息 的 反映 精 
度 。 具 体 公式 如 下 : 


NDVI = ( g2 - pl)/(p2 *pl) (2) 
EVI 2 G[ (p2 -pl)/(p2 * C1 xpl - C2 xp3) +L] (3) 


式 中 :p1,p2 ,p3 分 别 对 应 MODIS 的 第 1 波段 (620 ~ 
670 nm) ,第 2 波段 (841 ~ 876 nm) 和 第 3 波段 
(459 ~479 nm) 反 射 率 。 系 数 工 为 冠 层 背 景 调整 参 
数 ,C1 ,C2 为 气 溶胶 电阻 项 系数 ,G 为 增益 因子 UU 
2.2 TVDI 原理 

土壤 水 分 的 多 寡 影 响 着 地 表 温 度 的 高 低 。 当 土 
壤 水 分 充足 ,植被 的 蒸 散 作用 增强 ,使 得 地 表 温 度 下 
降 。 相 反 ,土壤 水 分 较 少 时 ,植被 蒸 散 作用 降低 ,地 
表 温 度 升 高 。 土 壤 湿度 状况 与 地 表 温 度 之 间 存 在 明 
显 的 负 相 关 性 。PRICE 等 ”发现 遥感 影像 所 表现 


出 的 温度 一 植被 指数 特征 空间 的 散 点 图 为 三 角形 。 
SANDBOIT 4&'" 对 简化 的 TANDT 三 角形 空间 进 
行 了 研究 ,提出 表示 土壤 湿度 的 温度 植被 干旱 指数 
TVDI 的 概念 ,定义 为 : 


TVDI = Ts A Ts VI) min (4) 
7 Ts( VI) = TA VI) min 

Ts( VI) max 一 Q max + Ds x VI (5 ) 

T, ( VI) min 一 d min + b min x VI (6 ) 


式 中 ;7, 代表 任意 一 点 像 元 的 地 表 温 度 值 , Ts 
(VI) 为 干 边 拟 合 方 程 ,Ts ( VI) 为 湿 边 拟 合 方 
程 ,a 和 5, 为 干 边线 性 拟 合 参 数 ,aji, 和 5, 为 湿 
边线 性 拟 合 参数 。 太 在 本 文中 表示 NDVI Ail EVI, 
2.3 改进 型 TVDI 的 构建 

2.3.1 mTVDI 模型 构建 过 程 TVDI 作为 一 种 水 
分 胁迫 指标 ,通过 LST 和 人 克之 间 建 立 T/VI 特征 空 
间 来 提取 干 湿 边 方程 ,进而 评估 研究 区 的 旱情 状况 。 
对 TVDI 模型 进行 改进 需要 考虑 数据 本 身 以 及 研究 
区 特有 的 地 理 环境 。 就 数据 源 本 身 而 言 ,遥感 数据 
源 会 受到 云 的 影响 使 得 模型 精度 大 大 降低 ,因此 需 
要 去 除 云 的 干扰 ;其 次 新 疆 干旱 区 面积 广大 ,地 势 起 
伏 较 大 ,地 表 温 度数 据 与 实际 地 表 温 度 之 间 易 产生 
偏差 ,也 会 降低 模型 精度 。 因 此 本 文 对 TVD 模型 
的 改进 主要 有 两 个 方面 :(1) 对 遥感 影像 进行 去 云 
处 理 。(2) 对 地 表 温 度数 据 进 行 地 形 校正 。 改 进 后 
的 TVDI 模型 即 为 mTVDI。 具 体 工作 包括 :(1) 利用 
MOD13A2 数据 提取 不 同时 相 NDVI 和 EV, 将 相同 
IAA LST 数据 和 NDVI,EVI 数据 按照 时 序 一 一 对 
应 。(2) 对 遥感 影像 进行 去 云 处 理 。(3 ) 对 LST 数 
据 进 行 地 形 校 正 。(4) 计算 T, 及 校正 后 T, 5 VI 
的 特征 空间 并 提取 干 湿 边 方程 。(5) 利用 不 同 干 湿 
边 方程 分 别 计 算 TVDI, , TVDI,, mTVDI, , mTVDI;。 
具体 构建 流程 如 图 2。 

2.3.2 去 云 处 理 ”由 于 气候 天 气 的 原因 ,大 部 分 
影像 在 获取 时 或 多 或 少 都 会 受到 云 的 影响 ,尤其 是 
大 范围 .长 时 序 光 谱 数 据 。 云 的 存在 使 得 遥感 影像 
在 某 些 区 域 存在 数据 缺失 ,大 气 背 景 条件 存 在 明显 
差异 ,同时 ,也 会 影响 到 地 表 温 度 的 反 演 精度 。 为 了 
保证 TVDI 模 型 反 演 精度 ,需要 对 影像 进行 去 云 处 
理 。 现 有 的 处 理 方式 有 两 种 :一 种 是 将 受 云 覆 盖 的 
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还 感 数据 获取 


| 
数据 预 处 理 


原始 TVDI 和 mTVDI 计 算 


图 2 mT VDI 模型 构建 


Fig.2 mTVDI model construction 


区 域 进行 数据 恢复 ,消除 云 的 影响 ; 另 一 种 是 检测 云 
的 存在 ,利用 掩 膜 进行 去 除 ,不 参与 模型 计算 。 为 了 
保证 反 演 精度 ,本 文 将 有 云 的 区 域 进行 掩 膜 计算 ,去 
除 云 对 TVDI 模型 精度 的 影响 。 

2.3.3 地 形 校正 利用 TYD 模 型 进行 旱情 监测 ， 
除了 受到 地 表 温 度 .植被 覆盖 度 的 影响 外 ,还 受到 地 
表 起 伏 状 况 ,土壤 质地 及 其 所 造成 的 大 气 辐 射 不 均 
衡 等 多 方面 的 影响 。 因 此 本 文 利 用 数字 高 程 数据 对 
地 表 温 度 进行 校正 ,以 消除 地 形 对 大 区 域 太阳 辐射 
差异 的 影响 。 公 式 如 下 : 


T, =T,+mxH 


(7) 


式 中 :有 代表 校正 后 的 地 表 温 度 T, 代表 原始 地 表 
温度 , 甩 为 高 程 值 ,m 为 高 程 修正 系数 。 

2.3.4 FRARR ， 干 湿 边 的 提取 是 模型 构建 的 
重要 内 容 。 干 湿 边 提取 即 利 用 植被 指数 数据 与 最 大 
最 小 地 表 温 度数 据 进行 线性 拟 合 ,获取 干 湿 边 方程 。 
本 文 提取 干 湿 边 的 方法 :V/ 值 范围 :0 ~1, 以 0.01+ 


由 表 2 可 以 看 出 ,改进 7;/V 干 湿 边 方程 拟 合 
效果 总 体 优 于 原始 Ts/ 和 干 湿 边 方程 拟 合 效果 。 
T,/EVI 和 改进 TsAEVI 拟 合 效果 优 于 T,/ NDVI 和 改 
进 T,/NDVI HA OR; M TEV 干 边 拟 合 方程 
R^ 在 0.8 以 上 , 湿 边 拟 合 方程 RR 在 0.5 以 上 , 拟 合 
效果 最 好 ,表现 出 了 一 定 的 优势 。 


3 结果 与 分 析 


3.1 模型 计算 与 分 析 
3.1.1 TVDI 计 算 结 果 对 比分 析 TVDI 的 取 值 范 
围 为 0 ~1, 当 TVDI ÉF 1 时 颜色 越 深 ,表示 旱情 愈 
加 严重 , 当 TVD ÉF 0 时 颜色 越 浅 ,表示 该 地 较为 
湿润 。 图 3 利用 TVDI,, TVDI, ,mTVDI,,mTVDI, 四 
种 模型 计算 出 的 2015 年 第 177 d 研究 区 旱情 空间 
分 布 情况 。 据 图 可 以 看 出 ,研究 区 TVDI 值 较 大 区 
域 主要 集中 在 盆地 地 区 ,盆地 地 形 使 得 大 气 不 易 散 
3A .蒸发 旺盛 .是 植 被 丰富 ,水 分 需求 量 巨 大 ,因此 干 
旱 程度 较为 严重 。7Y77 值 较 小 的 区 域 主要 分 布 在 
阿尔 泰山 脉 天 山 山脉 和 昆仑 山脉 及 其 周围 。 这 是 
由 于 这 些 区 域 在 这 个 时 期 冰雪 融 水 补给 较为 充足 ， 
所 以 表现 出 湿润 的 状态 。 四 种 模型 均 表 现 出 了 一 致 
的 旱情 分 布 特征 。 
3.1.2 TVDI 与 土壤 湿度 相关 性 比较 ”利用 TV- 
DI, ,TVDI, , mTVDI, , mTVDI, 模型 计算 结果 与 6 月 
份 土 壤 实 测 数 据 进 行 相 关 性 分 析 。 由 图 4 可 以 看 
H ,mTVDI, 与 实测 数据 相关 性 最 高 ,RE =0.74。 相 
比 其 它 模型 表现 出 了 一 定 的 优越 性 。 因 此 本 文选 取 
mTVDI, 作为 研究 区 旱情 监测 模型 ,对 2015 年 干旱 
区 植被 生长 季 旱 情 进行 监测 。 
3.2 mTVDI, 时 空 分 布 

图 5 为 2015 年 新 疆 植被 生长 季 mTVDI, 旱情 
监测 结果 。 从 空间 上 看 :(1) 旱 情 最 严重 的 地 方 主 
要 分 布 在 吐鲁番 盆地 及 哈密 盆地 西部 。 在 整个 生长 
季 当 中 , 叶 和 角 番 盆地 及 哈密 盆地 西部 都 处 于 极度 干 
旱 状态 。 主 要 原因 是 该 地 区 海拔 较 低 ,盆地 地 形 使 
得 热气 流 无 法 向 外 扩散 ,蒸发 强烈 ,旱情 最 为 严重 。 
(2) 旱 情 较为 严重 的 地 区 分 布 在 塔里木 盆地 ,该 地 


0. 02n(n = 1,2,…,99) 为 中 心 , 将 相同 VI 与 所 对 应 


区 沙漠 广 布 ,地表 薰 发 旺盛 , 且 四 周 有 山 脉 阻挡 ,水 


的 改进 后 的 Ts 的 最 大 值 和 最 小 值 对 应 ,最 后 采用 最 
小 二 乘法 拟 合 干 边 和 湿 边 。 利 用 R 语言 编程 ,自动 
提取 干 湿 边 。 植 被 生长 季 ( 第 97 ~ 305 d) 干 湿 边 方 
程 的 结果 如 表 2 所 示 。 


汽 难以 到 达 形 成 降水 ,致使 该 地 旱情 较为 严重 。 
(3) 旱 情 较为 严重 的 另 一 地 区 为 准噶尔 盆地 ,但 其 
旱情 状况 波动 较 大 。 在 第 129 d,161 d 和 177 d 旱 
情 较为 严重 ,而 其 他 时 间 旱 情 有 所 缓解 甚至 表现 出 


tió wy 


正常 的 状态 。 这 是 由 于 该 地 区 为 新 疆 重 要 的 农业 种 
植 区 ,人 为 灌溉 活动 较为 频繁 ,使 得 该 区 域 旱情 状况 
不 稳定 。(4) 旱 情 不 明 显 的 区 域 即 正常 状态 下 的 区 
域 主要 分 布 在 天 山 山脉 阿尔 泰山 山脉 和 昆仑 山脉 
附近 ,这些 地 方 在 植被 生长 季 时 期 既 有 雨水 补给 又 
有 高 山 冰 雪 融 水 补给 ,因此 并 不 会 出 现 明显 的 干旱 
状况 。(5 ) 旱情 不 明显 的 区 域 还 有 塔里木 河 周围 ， 
在 植被 生长 季 时 期 ,河流 对 周围 农田 和 树林 等 绿洲 
区 起 到 了 补给 灌溉 的 作用 。(6 ) 较为 湿润 的 区 域 分 
布 在 天 山 山脉 阿尔 泰山 山脉 和 昆仑 山脉 。 由 于 植 


相 比 较 可 以 看 出 :(1) 6,7,8 月 份 全 疆 降 水 偏 少 ,对 
应 的 mTVDI, 值 也 偏 高 ,全 疆 呈 现 出 严重 干旱 态势 。 
(2) 9,10 月 份 全 疆 降 水 偏 多 ,对 应 的 mTVDI, 值 偏 
低 ,全 疆 旱 情 有 所 缓解 。 

为 进一步 探究 新 疆 旱 情 与 降水 的 关系 ,本 文 对 
4 «10 月 mTVDI; 结果 与 降水 数据 进行 相关 性 分 析 。 
结果 如 图 7:4 月 尼 20.449 4,5 H R =0.415 9,6 
H R = 0.409 6,7 H R^ =0.431 4,8 H R =0. 43,9 
AR’ =0. 444 2,10 H R^ =0.4613。 由 此 可 以 看 
出 ,生长 季 各 月 TVDI 与 降水 之 间 存 在 负 相 关 关 系 ， 


被 生长 季 阶 段 温 度 较 高 ,冰雪 融 水 丰富 ,因此 表现 出 
较为 湿润 的 状态 。 
从 时 间 分 布 上 看 : (1) 第 161, 177, 193 ,209 , 
225 ,241 天 旱情 较为 严重 ,这 是 由 于 2015 年 伊朗 副 
高 正在 新 疆 上 空 稳定 停留 且 不 断 增强 所 致 。(2) 其 ”4 
他 时 间 段 新 疆 整体 相对 湿润 。 由 于 处 于 换季 时 期 ， 
中 ”研究 区 整体 降水 增多 ,m7VD1; 值 较 小 ,全 疆 早 情 
5 IRZ 
CO 3.3 TRMM 降水 数据 与 mTVDI, 相关 性 分 析 


其 结果 均 通 过 了 P <0. 01 显著 性 检验 。 综 上 ,新 疆 
旱情 与 降水 表现 出 一 定 的 相关 性 ,降水 的 多 少 影响 
着 新 疆 植被 生长 季 旱 情 的 广度 及 严重 程度 。 


疆 


结论 


为 了 提高 TVDI 模型 对 植被 生长 季 旱 情 监测 的 
敏感 性 ,提高 TVDI 模型 旱情 监测 精度 , 本文 利用 植 
被 指数 与 改进 后 的 地 表 温度 数据 构建 了 改进 型 的 早 
情 监 测 模型 mTVDI , 对 比分 析 TVDIL, , TVDI ; , mTV- 
DI, ,mTVDI, 四 个 模型 指示 旱情 的 优 劣 ,选取 最 佳 模 
型 对 2015 年 新 疆 干 旱 区 植被 生长 季 旱 情 进行 时 空 
分 析 并 利用 TRMM 降 水 数据 对 其 进行 验证 。 得 到 


co 图 6 利用 每 月 平均 每 小 时 降水 数据 乘 以 24 fi 
SS 算得 到 平均 每 月 每 天 降水 量 数据 。 根 据 各 月 平均 每 
天 降水 分 布 来 看 :6,7,8 月 降水 分 布 区 域 较 小 ,9 ,10 
月 份 降水 量 分 布 区 域 较 大 。 与 mTVDL, 时 空 分 布 图 


表 2 干 湿 边 方程 拟 合 方程 
© Tab.2 Dry and wet edge fitting equation 


一 原始 T,/NDVI 干 边 方程 R 湿 边 方程 R? 
C 第 97 天 y= -102. 389x +322. 43 0.813 7 y 254. 989x 4 223. 620 5 0. 463 8 
第 116 X y= -72.316x 4 332. 909 0. 894 7 y 2770. 537 2x +225. 261 3 0.633 9 
第 129 X y= -64.800 1x 4337. 700 2 0.923 8 y 268.915 2x +225. 878 4 0.8134 
第 145 X y = -62.147 8x +338. 749 4 0. 907 8 y 262. 822 1x +220. 172 2 0. 863 6 
第 161 天 y= -61.905x +344. 246 3 0. 854 2 y =51. 684 6x +225. 582 4 0. 808 5 
第 177 天 y= -49. 183x +339. 811 0. 863 1 y 251.178 6x +229. 148 4 0. 868 7 
第 193 X y= -51. 961 7x +342. 5534 0. 846 3 y =40. 735 7x +240. 95 0. 675 3 
第 200 天 y = -49. 987 9x +345. 148 3 0. 960 6 y 259.178 6x +234. 631 1 0. 855 
第 225 天 y= -40. 713 8x +336. 680 8 0.961 5 y 2773. 300 5x +230. 373 0. 849 3 
第 241 天 y= —38. 135 5x +330. 078 9 0.931 1 y 285. 971x + 227. 824 8 0. 776 
第 257 天 y= —51.897 2x +328. 308 3 0.758 5 y =75. 330 9x +230. 586 0.554 8 
第 273 天 y= -52.392 6x +320. 093 6 0. 801 8 y 272.284 9x +233. 311 1 0. 883 6 
第 289 天 y= -55.619 6x +311. 206 2 0. 761 9 y 229. 586 4x +241. 466 0.470 9 
第 305 X y = —90.498 6x +308. 063 7 0. 857 9 y 241. 236 4x +234. 459 2 0.577 2 
第 97 天 y= -67. 532 7x +322. 143 6 0. 845 4 y 231. 684x +230. 736 4 0. 501 9 
第 116 X y= -51.926 1x +332. 059 9 0.902 1 y 251. 342x +231. 032 9 0. 749 5 
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( 续 表 ) 

原始 Ts/EVI 干 边 方程 R 湿 边 方程 R? 
第 129 X y= —49.945 4x +337. 841 7 0. 898 7 y 243. 921x +233. 230 6 0. 737 1 
第 145 X y= -48.998x +338. 821 7 0. 889 y =51. 904x +226. 430 2 0. 895 
第 161 X y= —48. 621 4x +343. 684 8 0. 835 9 y 244. 334x +229. 115 8 0.8182 
第 177 天 y= —44. 289 7x € 341. 072 7 0. 877 1 y 242. ATI 3x +230. 954 4 0. 842 3 
第 193 天 y= —49. 176 4x +344. 857 5 0. 807 2 y =31.538 2x +240. 222 2 0.533 4 
第 209 X y= —46. 452 5x +347. 528 5 0. 909 2 y 245. 192x +234. 918 8 0. 795 4 
第 225 X y = —43. 085 8x +340. 691 9 0. 823 6 y 246.081 1x +234. 484 8 0. 898 8 
98 241 X y= -42.144 1x € 334. 2512 0. 795 8 y 243. 795 4x +234. 134 4 0. 734 1 
第 257 X y= —48. 654 4 «331.199 1 0. 802 6 y 239. 712 3x +235. 501 9 0. 743 9 
第 273 天 y= —40. 861 3x +320.981 1 0. 869 y =32. 771 2x +236. 918 3 0. 714 2 
第 289 天 y= -45.007 lx +315.936 8 0. 862 y =25. 699 7x +237. 21 0.582 2 
第 305 天 y = —40. 190 9x +303. 949 0. 864 y =23. 866 2x +235. 035 7 0.541 6 

改进 T,/NDVI 干 边 方程 Rm 湿 边 方程 Rm 
第 97 天 y= —68.395 2x +318.961 5 0. 837 4 y =27. 023x +259. 405 5 0.5434 
第 116 X y= -63.326 3x +331. 267 1 0. 880 4 y 221. 723 8x +264. 795 3 0.404 9 
第 129 X y= —-48. 732 1x +333. 657 0. 954 5 y =36. 784 3x +262. 161 6 0.855 8 
第 145 天 y= -51.031 2x +336. 224 0.9512 y 230. 026 8x +260. 674 7 0. 831 7 
第 161 天 y= —47. 429 8x +340. 361 5 0. 891 7 y 225. 780 8x +264. 100 9 0.779 3 
第 177 天 y= —39.578 2x +337. 462 1 0.927 1 y =28.548 lx +265. 721 7 0. 883 7 
第 193 X y= —42. 392 2x +339. 938 0. 932 y 220. 184 lx +275. 549 8 0. 609 4 
第 209 天 y = —45. 355 6x +344. 354 5 0.946 4 y =30. 756 6x +270. 765 9 0. 789 6 
第 225 天 y = -35.881 6x +335. 860 4 0.951 9 y 241. 071 2x +267. 526 5 0. 863 6 
第 241 天 y= —32.271 8x +329.047 1 0.928 8 y 242. 672 8x + 267. 673 7 0. 733 1 
第 257 天 y= -38.248 4x +325.321 8 0. 851 y 237. 585 5x +267. 5614 0.5315 
第 273 天 y= -47.003 8x +319. 31 0.903 1 y =38. 444x +260. 642 9 0. 865 5 
第 289 X y= -37.150 9x +310. 513 7 0. 8519 y 212.360 1x +271. 839 6 0. 434 1 
第 305 天 y= —52.514 7x +306. 694 3 0. 89 y 233.452 1x +259. 471 6 0.799 1 
改进 7s/EVI 干 边 方程 R 湿 边 方程 R? 
第 97 天 y= —45.028 6x +318. 463 1 0. 862 9 y =16. 761 4x +261. 889 0.676 4 
第 116 天 y= -42.915 5x +329. 723 1 0. 940 9 y =16. 843 1x +265. 364 3 0. 586 3 
第 129 天 y= —42. 304 3x +335. 525 0. 885 4 y =14.537 8x +267. 895 5 0.658 5 
第 145 X y= -42. 798 8x +337. 232 1 0. 908 1 y =22. 244 9x +263. 886 6 0.853 5 
第 161 天 y = —40. 736 8x +341. 42 0. 852 6 y =21.567x +265. 940 4 0. 852 6 
第 177 天 y= -38.307 2x +339. 618 2 0. 906 7 y 221. 667x +267. 122 8 0. 860 2 
第 193 天 y= —41. 691 9x +342.495 1 0.878 2 y =15. 127x +275. 271 8 0.578 4 
第 209 天 y= —42. 335 3x +346.511 8 0.895 7 y =26. 109 9x +269. 747 7 0. 833 6 
第 225 天 y= —37. 106 6x +338.991 0. 862 y =25. 366 4x + 269. 655 3 0. 880 9 
第 241 X y= —36. 148 8x € 332. 788 5 0. 836 9 y 219. 820 3x +271. 411 3 0. 709 4 
第 257 天 y= —37.859 4x € 328. 175 5 0. 901 5 y 217. 960 1x € 270. T75 2 0.643 8 
第 273 X y= —34. 499 9x +319. 424 1 0. 930 9 y =22. 418 6x +260. 981 6 0.857 9 
第 289 天 y = -33.404 9x+314.214 9 0. 966 y =12.524 lx +269. 129 6 0. 642 4 
第 305 天 y = —28. 788 4x +304. 976 6 0. 940 9 y 218. 467 7x +259. 854 0. 877 4 
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图 4 TVDI 与 实测 土壤 水 分 相关 性 比较 


Fig.4 Correlation between the TVDI and measured soil moisture 
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陈 丙 寅 等 ;基于 改进 型 TVDI 在 干旱 区 旱情 监测 中 的 应 用 研究 
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图 5 研究 区 2015 年 植被 生长 季 mTVDI, 时 空 分 布 
Fig.5 Temporal and spatial distribution of mTVDI, in the growing season of 2015 in the study area 


201908.00006v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


P 4 - 
tió wi 
z 80°E 90°E Zz = 80°E 90° E z 
3 l 名 a i N18 
Zz Z z 
©, 上 o Z o 
E + & 1g 
降水 量 /mm*d-! 
高 : 6.263 
1 000 km 
i 低 : 0.285 
80°E 90°E 80°E 90°E 
Z a E E E E E: T 
3 2 a 3 
Z z Zz 
o o zZ ae 
E E & S 
D 2 $ 2 
Z z z 
i: s S IE 
zZ zZ 
3 è 
图 6 研究 区 2015 年 4~10 各 月 平均 每 天 TRMM 

Z 
el = 降水 分 布 图 
+ 

3 Fig.6 Distribution of average daily TRMM precipitation 

mÓ 
eas in the study area from April to October in 2015 
0 1000 km = Nom 
———i| 
低 : 0.053 


80° 90°E 
以 下 结论 :(1) 经 过 地 形 校正 后 的 mTVDI, 对 干旱 “地 和 准 路 尔 盆地 为 两 个 干旱 中 心 , 旱 情 状 况 由 严重 
区 旱情 监测 效果 最 好 ,与 土壤 实测 数据 相关 性 最 高 ， “逐步 向 周围 山区 递减 至 湿润 状态 。(3 ) 从 时 间 上 
R'-0.74, (2) 新 疆 2015 年 旱情 分 布 以 塔里木 盆 看 ,6,7,8 月 份 全 疆 旱 情 较 为 严重 ,所 受 干 旱 影响 面 


陈 丙 实 等 :基于 改进 型 TVDI 在 干旱 区 旱情 监测 中 的 应 用 研究 
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积 较 大 。(4) 降水 与 新 疆 旱 情 具 有 一 定 的 相关 性 ， 
6,7,8 月 份 全 疆 降 水 偏 少 ,9 ,10 月 份 全 疆 降 水 偏 多 ， 
与 对 应 的 mTVDI, 结果 表现 出 一 致 的 变化 特征 ,由 
此 可 见 降水 多 窒 影 响 着 新 疆 干旱 区 旱情 状况 。(5) 
mTVDI, 与 TRMM 具有 一 定 的 相关 性 ,m7TVDI; 对 研 
究 区 旱情 状况 具有 一 定 的 指示 作用 。 

综 上 ,基于 增强 型 植被 指数 与 经 过 地 形 校 正 
后 的 地 表 温 度数 据 所 构建 的 m7TVDI; 模型 能 够 满 
足 干旱 区 植被 生长 季 旱 情 监 测 的 目的 。m7TVDI; 
模型 所 表现 出 的 旱情 状况 与 土壤 水 分 之 间 显 示 出 
较 强 的 相关 性 ,对 全 球 干 旱 区 旱情 研究 具有 很 强 
的 参考 意义 。 
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图 7 植被 生长 期 各 月 mTVDT 与 降水 量 相关 性 
Fig.7 Correlation between mTVDI, and monthly 
precipitation in growing season 
5 讨论 


mTTYD7 模 型 能 够 方便 高 效 的 反映 大 区 域 旱 情 
状况 ,但 还 存在 以 下 不 足 之 处 ,主要 表现 在 :(1) 般 
感 影像 质量 依然 是 影响 TVDI 模型 精度 的 重要 因 
素 。(2) 新 疆 地 区 面积 广大 ,土壤 质地 土地 利用 
类 型 的 不 同 所 表现 出 的 旱情 状况 也 不 尽 相 同 。 
此 接 下 来 的 工作 将 从 以 上 两 个 方面 进行 研究 。 
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Application of modified TVDI in drought monitoring in arid areas 


CHEN Bing-yin"^, YANG Liao', CHEN Xi’, WANG Wei-sheng! 
(1 Xinjiang Institute of Ecology and Geography ,State Key Laboratory of Desert and Oasis Ecology , Urumqi 830011 , Xinjiang , China; 
2 University of Chinese Academy of Sciences , Beijing 100049, China; 3 Research Center for Ecology and Environment of 
Central Asia , Chinese Academy of Sciences , Urumqi 830011 , Xinjiang , China) 


Abstract; The drought is one of the most influential natural disasters on a global scale. The soil desertification 
and salinization from drought can cause irreversible damage to the ecological environment. On basis of the prepro- 
cessing and cloud removal of MODIS data ,this paper modifies the TVDI model through terrain correction and con- 
structed an modified Temperature Vegetation Dryness Index ( mTVDI) for drought monitoring in arid regions of Xin- 
jiang , China ,and compared it with the traditional TVDI and verified it using soil measured data. The results showed 
as follows: (1) mTVDI , constructed by EVI ( Enhanced Vegetation Index) and corrected LST( Land Surface Temper- 
ature ) has the highest sensitivity to drought monitoring in arid area, and the correlation with measured soil moisture 
data is R^ :0. 74; (2) In 2015 ,the drought in Xinjiang was centered in the Tarim Basin and the Junggar Basin, and 
the situation was gradually decreased from the serious drought in the centers to a humid state in the mountainous ar- 
eas. From the time point of view ,the drought in Xinjiang was the most serious in June, July and August; (3) The 
spatial and temporal distribution characteristics of drought in Xinjiang based on mTVDI,, were compared and ana- 
lyzed with TRMM precipitation data. The result shows that the spatial and temporal distribution of drought was con- 
sistent between the two methods. The correlation between precipitation and mTVDI, in different time periods was 
high, and passed the P < 0.01 significance test. Above all, m7TVDI,, based on TVDI can effectively carry out drought 
monitoring in arid regions of Xinjiang with high precision. Therefore ,this will provide reference for the quantitative 
research on drought monitoring at large regional scale in the future. 


Key words: arid area; drought monitoring; TVDI; TRMM 


